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Abstract 

Effects of tlie Dirac sea on the excitation energy of the giant monopole states are 
investigated in an analytic way within the a — lo model. The excitation energy is 
determined by the relativistic Landau- Migdal parameters, Fq and F\. Their ana- 
lytic expressions are derived in the relativistic random phase approximation(RRPA) 
without the Dirac sea, with the Pauli blocking terms and with the full Dirac sea. 
It is shown that in the RRPA based on the mean field approximation the Pauli 
blocking terms should be included in the configuration space, according to the rela- 
tivistic Landau theory. In the renormalized RRPA, the incompressibility coefficient 
becomes negative, if NN excitations are neglected. 
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The relativistic mean field approximation(RMFA) neglects the Dirac sea in 
the description of the nuclear ground state. Recently, however, it has been nu- 
merically shown that in the relativistic random phase approximation(RRPA) 
built on the RMFA, the monopole states cannot be well described without 
the Pauli blocking terms which express transitions between the Dirac sea and 
the occupied Fermi sea [1]. If the blocking terms are neglected, the excitation 
energies of the monopole states in the RRPA are much lower than those in 
the time-dependent relativistic mean field approximation [2]. 

The purpose of the present paper is to show in an analytic way the role 
of the Dirac sea in the excitation energy of the monopole states. We will 
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discuss the monopole states of nuclear matter in terms of the Landau-Migdal 
parameters using the a — uj model. First, we will show that in the RRPA based 
on the RMFA, one should take into account the Pauli blocking terms in the 
configuration space. Second, the real effects of the Dirac sea will be discussed 
in the renormalized RRPA. It will be shown that NN states yield essential 
effects on the excitation energy through the Landau-Migdal parameter, Fq. 

The Landau-Migdal parameters, Fq and Fi, are obtained by the second deriva- 
tive of the total energy density with respect to the quasiparticle distribution. 
In the RMFA, they are given by [3,4] 

F-F F F- (^\ 



where we have defined 

F.^nJ^)\ F, = nJ^)\ (2) 
\msJ KrriyJ 

A^P = ^, v,^^, E,^ipl + M*r'. (3) 



In the above equations, gs and stand for the Yukawa coupling constants, mg 
and niy the masses of the a- and a;-meson, respectively, and pp and M* denote 
the Fermi momentum and the effective nucleon mass. Np and v-p represent the 
relativistic density of states at the Fermi surface and the relativistic Fermi 
velocity. The factor, a^, in Fq of Eq.(l) will play an essential role in later 
discussions, which is given by 



where 9p denotes the step function, 9{pf — \p\). 

In the relativistic model, the excitation energy of the monopole states is ex- 
pressed as [5], 

/ K V 



where ep denotes the Fermi energy and (r^) the root mean square radius of 
the nucleus. The incompressibihty coefficient, K, is expressed in terms of the 
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above relativistic Landau-Migdal parameters, 



Since pp is determined by the nucleon density, and ep is related to the nucleon 
binding energy, E^, and the free nucleon mass, M, 

eF = EB + M, (7) 



the excitation energy of the monopole state is a function of Fq and Fi. 

In order to see the effects of the Pauli blocking terms on the monopole states, 
we derive the Landau-Migdal parameters according to the RRPA. We calcu- 
late the longitudinal RRPA correlation functions with and without the Pauli 
blocking terms. By comparing them with the correlation function of the Lan- 
dau theory[6] , we will obtain the expressions of the Landau-Migdal parameters 
in each approximation. 

When following our previous papers [4,7], the mean field correlation function, 
IIh [4,7] , is given by the Fourier transform of the single-particle Green function, 

UH(A,B;k) = d^pTr[rAGH{p + k)rBGH{p)], (8) 



where k denotes the four-momentum, [ko, k), and A and B the Fourier trans- 
form of the external field expressed with the mean field, ipu{x), 

A{k) = J d^xex.p{ik • x)'iljH{x)rAil^H{x), (9) 



being some 4x4 matrices. The sum of the ring diagrams in the RRPA for 
the a — (J model is described as [4,7], 

6UnPA{AB;k) = -^^H(AA";^)(C/L)a^.^H(A^S;A;), (10) 
dettyL 



where the contraction should be carried out with respect to the superfix and 
suffix, a,b — —1, 0, and A_i and Aq are given by Eq.(9) with F^ = 1 and 70, 
respectively. Moreover, Xs ^-nd Xv represent, 

Y =— ^ y=— ^ — (11) 
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The explicit form of the 2x2 matrix, Ui,, in Eq.(lO) depends on whether or not 
the Pauh blocking terms are included in the mean field correlation functions 
as discussed below. 



The Green function, Gu, is given by the sum of those for a single-particle, 
hole and antinucleon, 

GH^G^(l-ep)+Gi,ep + G^. (12) 
It is rewritten as a sum of the density- dependent and the Feynman part [7], 

Gh = Gd + Gp, Gb = 9p{G\y — Gp), G-p = Gp + (7^- (13) 

Hence, IIh is composed of the four terms hke Gt)Gb: G^Gp, GpGp) and GpGp 
[7]. In the RRPA based on the RMFA, the GpGp term is neglected, which is 
divergent, while in the previous calculations [4], the G^Gp and GpGp, terms, 
which contain the Pauli blocking NN excitations like G^OpG^, have been kept. 
Then we have obtained 



L 



Xs(i - xvHv) xsXvnsv ^ 



(14) 



where the mean field correlation functions are defined as 

n, = nH(A-i,A-i;A:), n, = nH(Ao,Ao;A;), (15) 
= nH(A_i, Ao; k) = nH(Ao, A_i; k). 

The Landau prescription of the correlation functions is obtained at the limit 
A; ^ 0. In this limit we have 

= (27r)=^ArF{(l - vlMx) - aj, (16) 
n, = {27r)^Np^x), n,. = (27r)3iVF(l - v^Y^^^x), (17) 

where $(a;) stands for the Lindhard function with x = ko/{\k\vp) [4]. Using 
these equations, we obtain the generahzed dielectric function from the factor 
of Eq.(lO) as 



^ -det[/L = (1 + a,F,) 1 1 + ( Fv - - v^F^xA $(x) I . (18) 
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In the Landau theory, Eq.(18) should be written as [6] 

^ -detC/L = cjl + ( Fo + -J^xA . (19) 



XsXv I V 1 + 3^1 



By comparing Eq.(18) with (19), we obtain the Landau-Migdal parameters 
which are the same as in Eq.(l). 

Next we investigate the role of the Pauli blocking terms in the Landau-Migdal 
parameters. We calculate the mean field correlation functions neglecting the 
Pauli blocking terms and taking the only particle-hole states. The calculation 
of the mean field correlation functions is a little different from the one in 
taking the Pauli blocking terms, since the correlation functions are not Lorentz 
covariant, and the continuity equation is provided in a different way, 

k^Uu{A,A'^;k) = {[Ao{k),A^k)]), (20) 



where A'^ is given by replacing P^ with 7'' in Eq.(9), and the r.h.s. is related 
to the expectation value of the ground state as, 

( |[Ao(fc),At(fc')]| ) = S{k-k'){[Aoik),A^k')]). (21) 



In including the Pauli blocking terms, the r.h.s. of Eq.(20) is always vanished, 
but in neglecting them, it is not for Ai^2,3- Hence, the relationship between the 
correlation functions due to the time- and the longitudinal component of the 
a;-meson is written in the frame, k — (|fe|, 0, 0), as 

nH(Ai, Ai; k) = y^nH(Ao, Ao; k) - a,, (22) 

where the additional term, Oy, comes from the r.h.s. of Eq.(20), 

av = ([Ao(fe),A^t(fe)])/|fc|. (23) 



Because of this fact, the back fiow effects due to the longitudinal a;-meson 
exchange on (^IIrpa are not simply normalized as Xv, and Ui^ in this case 
depends on Oy, 

I Xs(l + OvXy - XvHv) XaXyHs. , 
Ul- \ I , I24j 

XaXyHsv Xv(l - XsHs) 
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where we have defined 



^^=(2^1^' ' ^v = Xv-avX^. (25) 



At the hmit, A; — > 0, XvO^v becomes to be 



Moreover, rig of the present case has not the term, Cg, in Eq.(16), while riy and 
rigv are the same as in Eq.(17). As a result the generahzed dielectric function 
in this case is given by 



XsXv 

= (l - f ^v) {l + (^v - ^.(1 - 4> - 41^.=) *(.)) . (27, 

Finally comparison of the above equation with Eq.(19) provides us with the 
Landau-Migdal parameters in neglecting the Pauli blocking terms, 

Fo = - (1 - vl)F,, Fi = -vlF,. (28) 



The difference between Eqs.(l) and (28) is very clear. Fq and Fi in Eq.(28) 
have no denominator. In order to obtain the correct expressions of Fq and Fi 
within the RMFA, thus we need to include the Pauli blocking terms in the 
configuration space of the RRPA. 

In the Landau prescription, the denominators in Eq.(l) come from the self- 
consistent derivative of the effective mass and the baryon current with respect 
to the quasi-particle distribution, [3,5]. As to the effective mass in the 
RMFA, 




we have 



dM* 



1 

V 



El. 



M* 



El 



V 



(30) 



6 



V being the nuclear volume. The coefficient of dM* /drij in the r.h.s. yields agFg 
in the denominator of Fq, as seen in Eq.(4). In the Green function formalism, 
the effective mass of the RMFA is written as 



Hence the coefficient of dM* /drij in Eq.(30) is given by 



«si"s = -tA^ f— y np,„u(A; = 0), (33) 



On the other hand, we have shown that 

1 (9s 
(27r)3 \m,, 



where Ilpauii represents the Pauli blocking terms in Eq.(8) for = = 1. 
Thus it is seen that in the Landau prescription of the RMFA, the derivative 
of Gd includes implicitly the Pauli blocking terms. Indeed, we can prove that 



If we integrate the r.h.s. over pq, the Pauli blocking terms only remain. 

The same discussion is possible for the denominator of Fi. The self-consistent 
derivative of the current, j, as to the quasi-particle distribution provides [3], 

^ - - -v'F ^ (35) 



The coefficient of the second term yields the denominator of Fi. Using the 
Green function, the current of the RMFA is written as. 

Hence, the coefficient of the second term in Eq.(35) is given by 

(37) 



3 ^ Vmv/ J (27r)4 I ' 
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We can show that the above derivative of Gd required in the RMFA is ex- 
pressed by using the Pauh blocking terms as 



— ^ = GD7i<^N + GnIiGd - Gb^^- (38) 



Let us explore the effects of the Pauh blocking terms in more detail. The Pauli 
blocking terms reduce always the contribution of the cr-meson to Fq through 
the factor as in Eq.(l), since is positive, as seen in Eq.(4). On the other hand, 
the contribution of the a;-meson to Fq is not affected. Therefore, the value of 
Fq becomes always smaller, and the incompressibility, K, is reduced according 
to Eq.(6), when the Pauli blocking terms are neglected. On the contrary, the 
absolute value of Fi in Eq.(28), which has no denominator, becomes always 
larger, compared with the correct one in Eq.(l), so that K is enhanced in 
neglecting the Pauli blocking terms. 

We calculate the values of Fq, Fi and K using the following parameters [8] as 
an example, 

M = 939, rris = 520, rriy = 783 (MeV), 

= 109.626, gl = 190.431, (39) 



which reproduce the nucleon binding energy, Eb — —15.75 MeV at pp — 1-30 
fm-^ In this case, we have M* = 0.541M, vf = 0.451 and = 9.07 x 10"^. 
These values provide us with 

Fo = 0.569, Fi = -1.151, X = 544 MeV (40) 



in taking the Pauli blocking terms, and 

Fo = -0.368, Fi = -1.866, X = 357 MeV (41) 



in neglecting the Pauh blocking terms. Thus both Landau-Migdal parame- 
ters are strongly affected by the Pauli blocking terms, and, in particular, Fq 
changes its sign. Consequently, the value of K is fairly reduced in neglecting 
the Pauli blocking terms. In ^°^Pb, the reduction amounts to about 2.7 MeV 
for the present parameters. This fact may be observed in ref. [1] by numerical 
calculations. 

Now real effects of the Dirac sea should be explored with the fully renormal- 
ized RRPA, where the G^Gf term is also calculated in the Hartree correlation 
functions. Then, the Pauli exclusion principle in both the Fermi sea and the 
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Dirac sea is correctly taken into account. Such calculations based on the renor- 
malized Hartree approximation(RHA) have been done by the present authors 
in ref. [9]. For complete discussions, we quote those results here. The Landau- 
Migdal parameters in this case are given by 

i^O = F., - «renFs, F, ^ - , (42) 

where we have used the abbreviations: 

F^^nJ^Y, «ren = — (43) 

\mo/ l + Ogi^s + aD 

Formally the above equation is similar to Eq.(l), but the mass of the a;-meson 
is replaced by the bare mass, mo, in F^ and the Dirac sea yields an additional 
effect, in ctren [9]- 

As to Fi, essentially there is no additional effects from the Dirac sea. The 
renormalized correlation function from the GpGp term due to the a;-meson 
exchange disappears at the limit /c — > and has no contribution to Fi. Re- 
placement of rriy by mo comes from the fact that the cu-meson propagator is 
written in terms of the bare mass[9]. The value of Fi, however, depends on 
those of the Yukawa couphng constants used in the RHA. In order to reproduce 
the nucleon binding energy and the Fermi momentum mentioned before[9] , the 
RHA requires = 66.117 and g"^ = 79.927. These values give M* = 0.7306M 
and mo = 691.171 MeV, so that we obtain Fi = -0.620. 

On the other hand, Fo is strongly affected by the Dirac sea through od in 
Eq.(43). In using vp — 0.3502 of the present RHA, its value is much larger 
than that of OgFs, 

Od = 0.405, ttsFs = 0.0296. (44) 

This Dirac sea effects reduce strongly the contribution of the cr-meson to Fo 
through ctren, ^ud we have 

Fo = 0.676. (45) 

If Od were neglected, then the value of Fo would be —1.56, which means K < 0. 
This fact reflects that N-degrees of freedom play an important role to stabilize 
the nucleus in the RHA. 

Finally we give two comments. First, since the restoring force of the giant 
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quadnipolc states comes mainly from the distortion of the kinetic energy den- 
sity, its excitation energy depends on Fi[5], 



Thus, the Pauli blocking terms affect also the excitation energy of the quadru- 
pole states in the RRPA based on th RMFA. Moreover, the Pauli blocking 
terms should be taken into account in the description of the center of mass 
motion, which requires the correct Fi [5]. This fact has been observed in 
arguments on the spurious state of RRPA by Dawson and Furnstahl [10]. In 
the same way the Pauli blocking terms are necessary for discussions of the 
nuclear current or magnetic moments [11]. The isovector dipole states depend 
on Fi and F[, which also require the Pauli blocking terms. The detail will be 
published in a separate paper. 

Second, we note that Eqs.(5), (6) and (46) arc formally the same as those 
in nonrelativistic models, except for ep in their denominators [5]. In nonrel- 
ativistic models, it is replaced by the nucleon mass, M. However, they are 
related to each other as Eq.(7) through the nucleon binding energy which is 
negligible compared with the nucleon mass. Thus the relativistic correction 
to the excitation energies of the monopole and quadrupole states is less than 
1% for E^ = —15.75 MeV, if the values of the Landau-Migdal parameters in 
relativistic models are the same as in nonrelativistic ones. 

In conclusion, in the relativistic random phase approximation(RRPA) based 
on the relativistic mean field approximation, the Pauli blocking terms should 
be taken into account for consistent descriptions of the Landau-Migdal pa- 
rameters. The real effects of the Dirac sea on the excitation energy of the 
monopole states are studied in the RRPA built on the renormalized Hartree 
approximation. Then the nucleon-antinucleon states affect strongly the exci- 
tation energy of the monopole states through the Landau-Migdal parameter, 
Fq. The incompressibility coefficient becomes negative, if antinucleon-degrees 
of freedom are neglected. The Landau-Migdal parameter, Fi, is not affected 
formally by the renormalization. 
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